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A structural model for intermediate sized silicon clusters is proposed that 
is able to generate unique structures without any dangling bonds. This struc- 
\f} \ tural model consists of bulk-like core of five atoms surrounded by fullerene-like 

surface. Reconstruction of the ideal fullerene geometry results in the forma- 
tion of crown atoms surrounded by 7r-bonded dimer pairs. This model yields 
unique structures for Si33 , Si3g , and Si45 clusters without any dangling bonds 
, and hence explains why these clusters are least reactive towards chemisorp- 

<^> \ tion of ammonia, methanol, ethylene, and water. This model is also consis- 

tent with the experimental finding that silicon clusters undergo a transition 
from prolate to spherical shapes at Si27- Finally, reagent specific chemisorp- 
tion reactivities observed experimentally is explained based on the electronic 
structures of the reagents. 
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Smalley and co-workers carried out extensive studies on chemisorption of various reagents 
on intermediate sized silicon clusters Q. These studies revealed that the reactivity rates 
for ammonia (NH 3 ) chemisorption vary by over three orders of magnitude as a function of 
cluster size, with 21, 25, 33, 39, and 45 atom clusters being particularly unreactive. Clusters 
containing more than forty seven atoms are highly reactive and they do not display strong 
oscillations in reactivities as a function of the cluster size. Similar results were obtained 
with methanol (CH3OH), ethylene (C2H4), and water (H2O). On the other hand, nitric 
oxide (NO) and oxygen (O2) were found to react readily with all clusters in this size range. 

Several structural models have been proposed to explain these experimental data 
However, these models do not yield unique and consistent structures for different cluster sizes 
and each cluster structure has to be obtained independently. Furthermore, these models do 
not satisfy the essential criterion that the structures of the unreactive clusters should not 
have any dangling bonds. 

In this Letter we propose a consistent model that generates the structures of the inter- 
mediate sized unreactive silicon clusters in a systematic way. The structures thus generated 
are unique for Si33, Si3g, and Si45 clusters. Furthermore, these clusters do not have any dan- 
gling bonds and hence explains why these clusters are unreactive for chemisorption. Finally, 
our model is consistent with the experimental finding of Jarrold and Constant that silicon 
clusters undergo a shape transition from prolate to spherical shapes at Si 27 0. 

Our structural model for silicon clusters consists of bulk-like core of five atoms surrounded 
by fullerene-like surface. The core atoms bind to the twelve surface atoms, thus rendering 
them bulk-like with four-fold coordination. The surface then relaxes from its ideal fullerene 
geometry to give rise to crown atoms and dimer pairs. The crown atoms are called adatoms 
in surface science literature || . These crown atoms are formally three-fold coordinated and 
possess one dangling bond each. The dimer pairs are also formally three-fold coordinated, 
but they eliminate their dangling bonds through local n bonding. This model yields struc- 
tures with seventeen four-fold coordinated atoms, four crown atoms, and the rest dimer 
pairs. The essential feature of this construction is that the bulk-like core and fullerene-like 
surface make these structures stable. This model is applicable to clusters containing more 
than twenty atoms. 

Unlike carbon, silicon does not form strong delocalized n bonds. Consequently, fullerene 
cage structures [p],|T0| , which require strong delocalized 7r bonds, are energetically unfavor- 
able for silicon. For this reason, silicon clusters prefer as few surface atoms as possible. 
Nonetheless, the fullerene geometry for the surface, consisting of interlocking pentagons and 
hexagons, gives special stability to the surface atoms. Furthermore, since delocalized 7r 
bonding is not favorable in silicon, we expect the surface atoms to relax from their ideal 
fullerene geometry to allow for dimer formation through strong local n bonding. Our model 
accounts for all these facts. 

The 5-atom core in our model has the exact structure of bulk silicon with one atom in 
the center bonded to four atoms arranged in a perfect tetrahedral symmetry. There are 
numerous ways to orient the 5-atom core inside the fullerene cage and bind it to the surface 



atoms. Furthermore, structural isomers may exist for fullerenes of any size |10|Jll| . Thus it 



is possible to use this model to generate all possible structural isomers for any odd numbered 
intermediate sized cluster. However, a particular orientation of the 5-atom core inside the 
fullerene cage will yield structures with maximum number of dimer pairs and least number 
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of dangling bonds. Such isomers are likely to be most stable. 

The 28- and 40-atom fullerenes are the only ones belonging to the symmetry group in 
the 20-60 atom size range |L(J . We generate the Si 33 and Si4 5 structures by inserting the 5- 
atom core inside the 28- and 40-atom fullerene, respectively. We orient the 5-atom pyramid 
in such a way that the central atom (A, violet), the apex atom (B, blue) of the pyramid, and 
the crown atom (C, red) on the surface lie on a line. The C atom is the central atom shared 
by three fused pentagons. This atom is surrounded by three other surface atoms (D, green), 
which relax inwards to bind to the four core atoms B. This relaxation motion is identical to 
that necessary to form the 2x1 reconstruction on the bulk Si(lll) surface |T2,r3]. With 



an activation barrier of only m 0.01 eV [1^,1^] and gain of 2.3 eV/bond due to a bond 
formation between B and D fl3]Jl5j , such a relaxation of fullerene surface is feasible even at 
100 K. Finally, the remaining surface atoms (E, orange) readjust to form as many dimers as 



possible. This relaxation is similar to that on Si(100) surface yielding dimer pairs [|T~3| , [T6 



The dimers are cr-bonded pair of atoms whose dangling bonds are saturated through strong 
local 7r bonds. Because of the tetrahedral symmetry of the core as well as of the reconstructed 
fullerene cage, the final Si 33 and Si45 structures also have tetrahedral symmetry. The Si39 
structure is also generated in this way, starting from the 34-atom fullerene cage and stuffing 
five atoms inside it. However, the structure thus generated has only C 3 „ symmetry, because 
of the lower symmetry of the fullerene cage [1C]. 

The Si 33 , Si 3g , and Si 45 cluster structures thus generated are displayed in Fig. 1. These 
structures have seventeen four-fold coordinated atoms, four crown atoms, and variable num- 
ber of dimer pairs. The Si 33 structure has six dimer pairs, whereas Si45 has twelve dimer 
pairs forming four hexagons. On the other hand, the surface of the Si 3 g cluster consists of 
nine dimer pairs, three of which form a hexagon. 

We constructed Si 35 and Si 43 clusters also using the proposed model, but found that these 
structures do not possess the characteristic crown-atom-dimer pattern found in unreactive 
clusters. Consequently, these structures possessed large number of dangling bonds. This 
explains the highly reactive nature of these clusters. 

In principle, the proposed model can be used to construct Si 2 5 structure also. However, 
Si 2 o cage is too small to accommodate five atoms inside it without undue strain. Conse- 
quently, Si 25 will prefer a different geometry. Indeed, there is experimental evidence that 
clusters smaller than Si 2 7 do not possess spherical shapes characteristic of larger clusters; 
instead they seem to prefer prolate shapes J5J . Our inability to generate a compact structure 
for Si 2 5 explains this experimental observation also. 

We verified our model by constructing structures of all the clusters discussed here using 
the molecular modelling kits |T7|. These structures unambiguously demonstrated that only 
33-, 39-, and 45-atom clusters possess the crown-atom-dimer pattern on their surfaces. We 
then calculated the atomic coordinates of these clusters from the tetrahedral geometry of the 



5-atom core and the known geometries of the fullerene structures |T0| , assuming that all bond 
lengths are equal to the bulk value of 2.35 A. The structures thus generated are displayed in 
Fig. 1. We then relaxed these structures by carrying out molecular dynamics at 100 K using 
the Kaxiras-Pandey potential |18| . The final structures obtained from these simulations are 



nearly identical to the initial ones. This indicates that the proposed structures are locally 
stable. 

The computer generated structures displayed in Fig. 1 revealed that the crown atoms 
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are able to form a fourth bond to the core atoms B, thus rendering the B atoms formally 
five-fold coordinated. The B-C bond arises from the back donation of the electrons from C 
to B and it weakens the neighboring bonds through electronic repulsion. The dangling bond 
on the crown atom is thus eliminated and these magic number clusters become unreactive. 



The classical potentials available at present [|Iq , |i9| are most appropriate for bulk silicon 
and related structures. Such potentials may not be suitable for describing unusual bonding 
patterns, such as the five-fold coordinated silicon atoms found in these clusters. Conse- 
quently, we cannot use the available classical potentials to prove that the proposed structures 
correspond to the lowest energy structures. Finally, structural determination based on the 
first principles electronic structure calculations are extremely demanding computationally at 
present for such large clusters, especially if adequate basis sets and grid sizes are employed 
and the calculations are converged to better than 0.01 eV accuracy. Furthermore, such 
complicated calculations do not provide the conceptual framework for understanding clus- 
ter properties. On the other hand, our simple physically motivated stuffed fullerene model 
yields insights into the nature of bonding in silicon clusters and explains the experimental 
trends in reactivities. 

There exist several competing structural models |2|-§[ of silicon clusters that attempt to 
explain the experimental reactivity data of Smalley and co-workers M. For example, Kaxiras 
has proposed structures of Si 33 and Si 45 clusters based on the reconstructed 7x7 and 2x1 
surfaces of bulk Si(lll), respectively Q. However, this model does not explain the reactivity 
data since bulk surfaces are highly reactive. Furthermore, it is inconsistent that the surface 
of Si45 should be the metastable 2x1 surface rather than the more stable 7x7 surface 



T3| . Finally, the Si45 structure of Kaxiras has forty dangling bonds, which make it highly 



reactive, contrary to the experiments. To overcome this discrepancy between experiment 
and theory, Kaxiras has proposed that the dangling bonds on Si 45 form 7r-bonded chains, 



analogous to the 2x1 reconstruction of the bulk Si(lll) surface | 20| . However, silicon 
favors strong local n bonds over delocalized 7r-bonded chains. This is the reason why the 2 
x 1 reconstruction, involving 7r-bonded chains p0], is metastable with respect to the locally 



7r-bonded 7x7 reconstruction on the bulk Si(lll) surface fl3fl . For the same reason, the Si 45 
structure proposed by Kaxiras is a metastable and highly reactive structure. Indeed, Jelski 
and co-workers disputed the Kaxiras model of Si45 by constructing alternative structures for 
Si45 that are lower in energy, but do not possess any of the features of the reconstructed 
bulk surfaces [|]. 

The bonding characteristics of silicon differ in subtle ways from that of carbon. In 
carbon, delocalized 7r bonds are favored over local n bonds, whereas the opposite is true in 
silicon. For this reason, graphite is the most stable form of carbon at room temperature 
and atmospheric pressure, but not the graphite form of silicon. Likewise, the bulk (111) 
surface of the diamond form of carbon exhibits the 2x1 reconstruction, but not the 7 
x 7 reconstruction [pTp2f . These examples, illustrate how subtle differences in bonding 
characteristics determine possible crystal structures and surface reconstructions. The same is 
true of clusters and the models of cluster structures should account for these characteristics. 
Our model of silicon clusters accounts for these facts by focussing on structures that are able 
to form maximum number of a bonds and eliminate their surface dangling bonds through 
local 7r bonding. 

Our structure for Si 33 is identical to that proposed by Kaxiras H and Patterson and 
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Messmer . This structure has been shown to be locally stable J23| . But our Si45 structure 



is different from that of Kaxiras Q . However, we can generate the Si45 structure of Kaxiras 
by stuffing one atom inside a 44-atom fullerene cage and allowing for the reconstruction of 
the fullerene surface. Thus our model is very general, subsuming the Kaxiras model within 
it. 

The reactivity patterns of NO and O2 are different from those found for NH3, CH3OH, 
C2H4, and H 2 [[!]]. This may be explained based on the ground state electronic structures 
of these reagents. NH 3 , CH 3 OH, C2H4, and H 2 in their ground states have closed shell 
electronic structure with all electrons paired. On the other hand, NO and 2 in their 
ground states are 2 U g and 3 S~, possessing one and two unpaired electrons, respectively 
24fl . Consequently, NH3, CH3OH, C2H4, and H2O can chemisorb only at those sites where 



excess electron density is present due to dangling bonds. Such a selectivity gives rise to 
highly oscillatory pattern in the reactivities, because the number of dangling bonds varies 
as a function of cluster size. The magic number clusters are unreactive because they do not 
possess any dangling bonds. On the other hand, NO ( 2 II 5 ) and 2 ( 3 S~) can chemisorb 
anywhere, because these reagents carry the necessary dangling bonds for instigating the 
reaction anywhere on the cluster surface. Hence, NO and 2 readily react with all clusters 
and do not display the oscillatory pattern in their chemical reactivities. This explains the 
reagent specific chemisorption reactivities observed experimentally [|l|]. 

The magic number clusters are not completely inert towards the closed shell reagents 
0]. These clusters are more reactive towards NH3, CH3OH, and H 2 than towards C 2 H4. 
This is because NH3, CH3OH, and H 2 have lone pairs on either nitrogen or oxygen and 
these lone pairs have a small probability of instigating reaction on the cluster surface. A 
lone pair is a pair of electrons that is not part of a bond. C 2 H 4 does not have any lone 
pairs and hence the magic number clusters are quite unreactive towards this molecule. The 
electronic the structure of reagents thus explains even subtle variations in the reactivities of 
magic number clusters towards a group of related reagents. 

In summary, we propose a structural model for the unreactive silicon clusters containing 
more than twenty atoms. This model consists of bulk-like core of five atoms surrounded by 
reconstructed fullerene surface. The resulting structures for Si 33 , Si 39 , and Si 45 are unique, 
have maximum number of four-fold coordinated atoms, minimum number of surface atoms, 
and zero dangling bonds. Such unique structures cannot be built for other intermediate sized 
clusters and hence they will have larger number of dangling bonds. This explains why Si3 3 , 
Si39, and Si45 clusters are least reactive towards closed shell reagents ammonia, methanol, 
ethylene, and water [|TJ. Our model also indicates that Si 25 cluster cannot be formed in a 
spherical shape. This result is consistent with the experimental finding that silicon clusters 
undergo a shape transition from prolate to spherical shapes at Si 2 7 [0] . Finally, two distinct 
patterns of chemisorption reactivities observed experimentally are explained based on the 
electronic structures of the reagents. The reactivities of closed shell reagents depend on the 
available number of dangling bond sites, whereas the reactivities of free radical reagents are 
not so dependent. Consequently, only the closed shell reagents are sensitive to the cluster 
structure and hence exhibit the highly oscillatory pattern in reactivities as a function of the 
cluster size. 
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FIGURES 



FIG. 1. Structures of S133, S139, and S145 clusters obtained using the proposed stuffed fullerene 
model. The atoms in different chemical environments are colored differently from violet to red, 
while only the representative atoms are labelled from A to E. These clusters do not possess any 
dangling bonds and hence are least reactive towards the closed shell reagents ammonia, methanol, 
ethylene, and water. 
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